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Abstract 

We construct two-electron scattering states and verify their tensor product structure in the 
infrared-regular massless Nelson model. The proof follows the lines of Haag-Ruelle scattering 
theory: Scattering state approximants are defined with the help of two time-dependent renor- 
malized creation operators of the electrons acting on the vacuum. They depend on ground state 
wave functions of the (single-electron) fiber Hamiltonians with infrared cut-off. Convergence 
of these approximants as t — > oo is shown with the help of Cook's method combined with a 
non-stationary phase argument. Removal of the infrared cut-off in the limit t — > oo requires 
sharp estimates on the derivatives of these ground state wave functions w.r.t. electron and pho- 
ton momenta, with mild dependence on the infrared cut-off. These key estimates, which carry 
information about the localization of electrons in space, are obtained in a companion paper 
with the help of iterative analytic perturbation theory. Our results hold in the weak coupling 
regime. 

1 Introduction and results 

The last two decades have witnessed substantial progress in the mathematical understanding of 
scattering of light and matter in the framework of non-relativistic QED. These advances pro- 
vided a rigorous foundation for the physical description of the Compton [CFP07, CFP09, Pi03, 
Pi05, FGS04] and Rayleigh scattering [DG99, DG04, FGS02], involving one massive particle 
(the 'electron') and many massless excitations ('photons'). However, the case of Coulomb 
scattering, i.e., collisions of two electrons in the presence of photons, has remained outside 
of the scope of these investigations. This is a serious gap in our understanding, given the 
tremendous importance of this scattering process, ranging from Rutherford's discovery of the 
structure of atoms, to modern high-energy physics experiments. This paper is a first part of a 
larger investigation, whose goal is to put this important process on rigorous grounds. 

A general framework of scattering theory for several electrons in models of non-relativistic 
QED, which has its roots in Haag-Ruelle scattering theory [Ha58, Ru62], was known to ex- 
perts already more than three decades ago [Fr, A173]. However, a complete construction of 
scattering states was only possible in the presence of a fixed infrared cut-off (or non-zero pho- 
ton mass), due to severe infrared- and infraparticle problems. These difficulties were overcome 
at the single-electron level with the help of a novel multiscale technique [Pi03, Pi05], but the 
case of several electrons remained open to date. In this paper we provide a construction of 
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two-electron scattering states in the massless, infrared-regular Nelson model. Although the 
infraparticle problem does not arise in this situation, the infrared structure of the model is non- 
trivial and requires major refinements of the multiscale technique. To clarify the origin of these 
new difficulties and describe our methods to tackle them, let us now explain in non-technical 
terms the main steps of our analysis. 

Let H be the Hamiltonian of the translationally invariant, infrared-regular, massless Nelson 
model stated in (1.24) below. It describes second-quantized non-relativistic particles, which 
we call electrons though they obey the Bose statistics, interacting with massless scalar bosons, 
which we call photons. This Hamiltonian is a self-adjoint operator on the physical Hilbert 
space *H which is the tensor product of the electron Fock space r(f) e ) and the photon Fock 
space r(f)f). As H preserves the number of electrons, we can restrict it to the one-electron 
subspace obtaining the Hamiltonian H^ l \ This Hamiltonian has the standard decomposition 
into the fiber Hamiltonians at fixed total momentum H p , which are operators on the Fock 
space. In the infrared-regular case (1/2 > a > in (1.21) below), the fiber Hamiltonains 
have (normalized) ground states ijjp, (for P in some ball S centered at zero), corresponding to 
eigenvalues Ep. Superpositions of such ground states of the form 

th -irj d 3 Ph(P)if, P , heC 2 (S), (1.1) 

give physical single-electron states in *K. (Here n* is the standard identification between the 
fiber picture and the physical picture). We note that the time evolution of if/}, is given by 

e- iHt ^ h = \j/ ht , (1.2) 

where h t {P) := e~ lEpt h{P). In heuristic terms, to construct a scattering state of two electrons 
one has to 'multiply' two single-electron states in such a way that the result is a vector in *K. 
Haag-Ruelle scattering theory is, in essence, a prescription to perform such a multiplication. 
In the context of the Nelson model one could attempt to implement this construction in the 
following way: Let {f™ } me N t> e tne m-photon components of if/ p. Then (1.1) can be rewritten 
as follows 

oo 

if/ h = V — d 3 p d 3m k h(p)f™(h k m )a\h) . . . a*(k m ) r}*(p - k)Cl, (1.3) 
to 

where Q. e *H is the vacuum vector, rf{p), a*{k) are the electron and photon creation operators 
and k := k\ + ■ ■ ■ + k m . This suggests the following definition of the renormalized creation 
operator 

°° 1 r 

fr(h):=V— d 3 pd 3m kKp)C(k l ,...,k m )a*(k l )...a*(k m )r 1 *(p-k), (1.4) 
to Vm!J 

which creates the physical single-electron state i/rj, from the vacuum 1 . Now, given /ji,/j2 £ 
Cq(5) with disjoint (velocity) supports, the scattering state h describing an asymptotic 
configuration of two independent electrons ^ , \ph 2 is given by 

% h = lim e iHt f,*(h u )f,*(h 2 ,)Q, (1.5) 

l We found this definiton of the renormalized creation operator and a construction of scattering states of several 
electrons in the Nelson model with a fixed infrared cut-off in unpublished notes of J. Frohlich [Fr]. Such renormalized 
creation operators appear also in a work of S. Albeverio [A173]. 
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if the approximating vectors on the r.h.s. are well defined and the limit exists. 

Unfortunately, regularity properties of the m-photon components f™ ofijjp, which control 
the rate of convergence in (1.5), are difficult to obtain due to the fact that the eigenvalue Ep 
is located at the bottom of the continuous spectrum of H p . Therefore, we will not construct 
scattering states with the help of formula (1.5), but instead we will use more tractable approx- 
imating sequences: Let Iry be the fiber Hamiltonians with the infrared cut-off cr, defined 
precisely in (1.30). Their ground states \j/p j(r , corresponding to eigenvalues Ep^, satisfy 

lim^ Pi0 - = tf/p. (1.6) 

(T-»0 

Let {/p^JmeNo be the m-photon components of tAp.o-. The renormalized creation operators with 
the infrared cut-off cr have the form 

oo 

fc(fc):=Y-= d 3 pd 3m kh(p)f™{k u ...,k m )a\k0...a*{k m )7 ] *{p-k). (1.7) 
We will use them to construct the scattering states as follows 

Km ■■= £m e^e ( (fci,*:e t (M«> a-») 

where cr, = tc/t r , y > 4 and k > is the ultraviolet cut-off. The existence of this limit and its 
properties, which allow to interpret it as an asymptotic configuration of two electrons, consti- 
tute the main result of this paper (Theorem 1.3 below). The main steps of the proof, which rely 
on Cook's method combined with non-stationary phase analysis, are given in Section 2. We 
conjecture that the limit (1.5) also exists and coincides with (1.8). Supporting evidence for this 
conjecture comes from algebraic quantum field theory, where the existence of scattering states 
in the infrared-regular situation has been proven without introducing cut-offs [Dy05]. We do 
not expect any problems with generalizing our result to an arbitrary number of electrons or 
with changing the statistics of our 'electrons' from Bose to Fermi. However, our aim here is 
not to cover the most general situation, but rather to prepare grounds for our future investiga- 
tion concerning Coulomb scattering in the infrared-singular case. The additional complications 
coming from the presence of infrared photon clouds should first be tackled in the case of two 
electrons. 

Let us conclude this introductory discussion with some more technical remarks. We point 
out that the approximating sequence (1.8) has important advantages over (1.5): Due to the fact 
that Ep^ is an isolated eigenvalue, the relevant regularity properties of the eigenvectors i/r Pcr 
and their m-photon components f™ can be obtained using the analytic perturbation theory. 
Their behaviour in the limit cr — » can be studied with the help of the iterative multiscale 
technique. As this analysis is presented in a separate paper [DP12], it suffices to indicate here 
the new spectral difficulties encountered at the level of two-electron scattering. 

We recall that already in [Fr73, Fr74] a formula for p* was derived, which in the case of 
m = 1 has the form 

where the form-factor vZ. is given by (1.31). It is not difficult to obtain from this formula that 

cv^Qc) 

|/^)| <^-, (1.10) 
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where the constant c is independent of cr. However, to prove convergence of (1.8) via the 
non-stationary phase method, we also need bounds on derivatives dpi/ft , dpidpjfp^ with mild 
dependence on cr. To indicate various difficulties that have to be tackled to obtain such bounds, 
let us differentiate (1.9) w.r.t. P: 

where Pf is the photon momentum operator. As for the second term on the r.h.s. of (1.11) the 
main difficulty is to find an estimate on the vector dpuj/p^ with suitable dependence on cr. We 
recall that the existing bounds from [Pi03] give only Holder continuity of P i-> i//p tCr , uniformly 
in cr. By a refinement of the argument from [Pi03], we obtain in [DP12] the following bound, 
stated precisely in Proposition 1 . 1 below, 

Wdpn/fpJl < (1.12) 

cr A Q 

where c is independent of cr, A® > is the maximal admissible value of the coupling constant 
and the function Aq i-» 0% satisfies lim^^o 8^ - 0. Although the bound is not uniform in cr, 
such mild dependence on the infrared cut-off is sufficient for our purposes. Denoting by // the 
second term on the r.h.s. of (1.11), and proceeding as in the derivation of (1.10), we obtain 

1 cv£(Jfc) 

|//|<-i-^— . (1.13) 

As for the first term on the r.h.s. of ( 1. 1 1), which we denote by /, a crude estimate gives 

cv1.{k) 1 cvl(k) 

|/|<-^<-^A (1.14) 
\k\ 2 cr \k\ 

due to the presence of the additional resolvent and the support properties of v'L. This bound 
does not suffice for the purpose of constructing scattering states. With the help of the multiscale 
analysis we improve it to 

1 cvZ(k) 

\I\<^-^—, (1.15) 



which is only slightly worse than (1.10). Altogether we get 

1 cvZ(k) 



\d P .f P(T (k)\<— a — (1.16) 

■ u Ar 1 ' ' 1 



and an analogous bound for the second derivative. These bounds, and their counterparts for 
fp ,m> 1, constitute the main result of [DP 12], stated also in Theorem 1.2 below. 

This paper is organized as follows: In Subsection 1. 1 we recall the definition of the Nelson 
model with many electrons. In Subsection 1.2 we state the relevant results concerning spectral 
theory, which are proven in a separate paper [DP12]. In Subsection 1.3 we state the main 
result of the present paper which concern scattering theory of two electrons in the infrared- 
regular massless Nelson model. Section 2 presents the main steps of the proof and in the 
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subsequent sections we provide the necessary ingredients. Section 3 is devoted to vacuum 
expectation values of the renormalized creation operators. In Section 4 decay properties of 
these vacuum expectation values are derived with the help of the method of non-stationary 
phase. Important input in this section are our spectral results proven in [DP12] and summarized 
in Subsection 1.2. The more technical part of our discussion is postponed to appendices. 
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1.1 The model 

We consider an interacting system of massive spinless bosons, which we will call 'electrons', 
and massless spinless bosons, which we will call 'photons'. Let f) e := L 2 (M 3 ,d 3 p) and f)f := 
L 2 (R 3 ,d 3 k) be the single-electron and single-photon spaces, respectively, and let T(\) e ) and 
r(bf) be the corresponding symmetric Fock spaces. The (improper) creation and annihilation 
operators on r(f) e ) (resp. r(t)f)) will be denoted by rf{p),rj(p) (resp. a*(k),a(k)). They satisfy 
the canonical commutation relations: 

Vr}(p),rf(p')} = S(p-p'), [T](p),v(p')] = tfipWip')] = 0, (1-17) 
[a{k),a{k')] = 6{k - k'), [a{k),a{k')] = [a{k),a{k')] = 0. (1.18) 

The free Hamiltonians of electrons and photons are given by 

H e := J d 3 P n(p)r]*(p)T](p), Hf.= J d 3 kaj(k)a*(k)a(k), (1.19) 



where £l(p) = \ and co(k) = \k\. We recall that these operators are essentially self-adjoint 
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on C e , Cf, respectively, where C e /f c r(b e /f) are dense subspaces consisting of finite linear 
combinations of symmetrized tensor products of elements of C ( ^°(R 3 ). 

The physical Hilbert space of our system is *K := r(l) e ) ® T(bf) and we will follow the 
standard convention to denote operators of the form A ® 1 and 1 ® B by A and B, respectively. 
The Hamiltonian describing the free evoluton of the composite system of electrons and photons 
is given by 

H ir :=H e + H f (1.20) 

and it is essentially self-adjoint on C := C e ® Cf. Now let us introduce the interaction between 
electrons and photons. Let A > be the coupling constant, k — 1 be the ultraviolet cut-off 2 
and let 1/2 > a > be a parameter which controls the infrared behavior of the system. Given 
these parameters, we define the form-factor 

Va(k) := A r , (1.21) 

(2|*|)5 



2 We set k — 1 to simplify the proofs of Proposition 1.1 and Theorem 1.2, given in the companion paper [DPI 2]. In 
the present paper we will write k explicitly. 
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where %k e C ( ^°(R 3 ) is rotationally invariant, non-increasing in the radial direction, supported 
in S K and equal to one on B(i- So ) K , for some fixed < eo < 1- (We denote by S r the open ball 
of radius r centered at zero). 

The interaction Hamiltonian, defined as a symmetric operator on C, is given by the follow- 
ing formula 



d 3 pd 3 kva(k)j]*(p + k)a(k)Tj(p) + h.c. (1.22) 



For future reference we denote by be the first term on the r.h.s. of (1.22) and set H? = (H*)*. 

As indicated in [Fr73, Fr74], the full Hamiltonian H = Hf T + Hi can be defined as a self- 
adjoint operator on a dense domain in ( H. For the reader's convenience we outline briefly 
this construction: First, we note that both and Hi preserve the number of electrons. Let 
us therefore define *H^ n ' := T^ n '(l) e ) ® r(Bf), where T^"\l) e ) is the «-particle subspace of r(f) e ) 
and let H£ and H^ be the restrictions of the respective operators to "H^, defined of C in) :- 
C n t H (n) . As shown in Lemma A.l, using the Kato-Rellich theorem, each H {n) = H^ + H^ } 
can be defined as a bounded from below, self-adjoint operator on the domain of H b , which is 
essentially self-adjoint on C^ n \ Then we can define 

H-=Q)H M (1.23) 



as an operator on C. Since ± i have dense ranges on C^ n \ H ± i has a dense range on C, 
thus it is essentially self-adjoint on this domain. We stress that the above construction is valid 
both in the infrared-regular case (1/2 > a > 0) and in the infrared singular sitiation (a = 0). 
On C we have the following formula for H 

H= f d 3 pQ,(p)i]*(p)r](p) + J d 3 kco(k)a*(k)a(k) 

+ (J d 3 pd 3 kv^k)r]*(p + k)a(k)j](p) + h.c.). (1.24) 
It reduces to a more familial - expression on C (n) 

H („) = Y^)-JlL + J d 3 koj{k)a{k)a{k) + Y_ l J d 3 kva(k)(e ikx 'a(k) + e~ ikXi a(k)), (1.25) 



where Xj are the position operators of the electrons. Finally, we introduce the electron and 
photon momentum operators 

p' e - J d 3 pp i 7 1 *{p)7]{p), P\;= J d 3 kk i a\k)a{k), ie {1,2,3}, (1.26) 

which are essentially self-adjoint on C. We recall that H is translationally invariant, that is it 
commutes with the total momentum operators P', given by 

P i :=P i e + P i f , ie {1,2,3}, (1.27) 

which are also essentially self-adjoint on C. 
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1.2 Spectral theory 

In this section we collect some necessary information about the spectrum of the single-electron 
Hamiltonian H (l \ We recall that the analysis of the spectrum of H p was initiated in [Fr73, 
Fr74] and advanced in [Pi03] with the help of the iterative analytic perturbation theory. Further 
developments along these lines, which we describe in detail below, can be found in [DP12]. 
Interesting results on the spectrum of the Nelson model with a slightly different form factor 
were also obtained in [AH12] by different methods. 

We recall that due to the translational invariance can be decomposed into a direct 
integral of the fiber Hamiltonians H p as follows 

#(D _ n * J~® tfptfXl n . (1.28) 

Here II : — > L 2 (R 3 ; r(bfi)) is a unitary map given by II := Fe' PfX , where x 1 are the 
electron's position operators and F is the Fourier transform in the electron's variables. The 
Hamiltonians H p ' are self-adjoint operators on the Fock space r(hfi), where bfi = L?(S?,d 3 k) 
is the single-photon space in the fiber picture. Denoting by b*(k) and b(k) the creation and 
annihilation operators on r(hfi) one easily obtains from (1.25) that 

H p 1} = -(P - P f f + H f + J d 3 kv w (k) (b(k) + b\k)). (1.29) 

As a tool to study the spectrum of this Hamiltonian, we introduce auxiliary fiber Hamiltonians 
with infrared cut-offs 

H pI = \( p ~ p f) 2 + H t + J d 3 kv^(k) (b(k) + b*(k)). (1.30) 
The form-factor v- is defined as follows 

a 

(21*1)2 

where < cr < K,X[o-,K)(k) := lsi-(*)Af*:(*)» i s the complement of the ball of radius cr and 1 A 
is the characteristic function of a set A. By replacing v„ with v'L in (1.22), we obtain the inter- 
action Hamiltonian //i j(T and the corresponding full Hamiltonian H a with an infrared cut-off. 
The restriction of H a to the single-electron subspace, denoted Hy, has a fiber decomposition 
into the Hamiltonians H), . 

P,(T 

Since we are particularly interested in the bottom of the spectrum of Hy and Hp^, let us 
define 

E P := inf cr(Hf), E PsfT := inf aiH^. (1.32) 

As the model is non-relativistic, we restrict attention to small values of the total momentum P 
at which the electron moves slower than the photons. More precisely, we consider P from the 
set 

5 :={PeR 3 ||P|<iW} (1.33) 

for some P ma x > 0. Since we work in the weak coupling regime, we fix some sufficiently 
small Aq > 0, and restrict attention to A € (0, Aq]. The parameters P max and Aq are specified 
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in Proposition 1.1 and Theorem 1.2 below. P max remains fixed in the course of our analysis. 
The maximal coupling constant Aq is readjusted only in the last step of our investigation - in 
Theorem 2.1 - to a new value which is denoted A' Q . 

In the following proposition we collect the results concerning Ep,Epo- which will be 
needed in the present investigation. Many of these properties are known, but several are new, 
as we explain below. 

Proposition 1.1. Fix < a < 1/2 and let P max = 1/6. Then there exists Aq > s.t. for all 
PeS-~ Sp max , A € (0, A ] there holds: 

(a) SbPh Ep is twice continuously differentiate and strictly convex. S 3 P \- > Ep t(r is 
analytic and strictly convex, uniformly in cr e (0, k]. Moreover, 

\Ep-E Pi(r \<co-, (1.34) 
\8?pEpA<c, \8PpE P ^\<c, \8$E Pja \<clo*»* (1.35) 

for multiindices f$j s.t. \f$j\ — j, je {1,2, 3}. 

(b) For cr > 0, Ep a- is a simple eigenvalue corresponding to a normalized eigenvector typ^, 
whose phase is specified in [DP 12]. There holds 

WtifypA ± c l^ i0 (!- 36 ) 

for multiindices f3 s.t. < |/?| < 2. 

(c) For a > 0, Ep is a simple eigenvalue corresponding to a normalized eigenvector ij/p. 
Moreover, for a suitable choice of the phase oftpp, 

Up-^pA\<cct". (1.37) 

The constant c above is independent of a, P, A, a within the assumed restrictions. Clearly, all 
statements above remain true after replacing Aq by some Aq e (0, Aq]. The resulting function 
Aq h-> <5t can be chosen positive and s.t. lim^ _,o^o = ^" 

In the light of Proposition 1.1 we can define the subspace of renormalized single-electron 
states 

<H 1(r : = {IT I d 3 Ph(P)if/p (T \he L 2 (R 3 , d 3 P),supph c S }. (1.38) 

In the case of 1/2 > a > we also set 'Hi := < H\, a -= l 0- 

Large part of Proposition 1 . 1 has already been established in the Nelson model or in similar 
models: The fact that 5 3 P \-> Ep is twice continuously differentiable and convex has been 
shown in non-relativistic and semi-relativistic QED in [FP10, KM12, BCFS07] and in the 
Nelson model with a slightly different form-factor in [AH 12]. The present case is covered by 
our analysis in [DP12]. The bound (1.34) can be extracted from [Pi03]. First statement in (b) 
has been established already in [Fr74]. Part (c) is implicit in [Pi03] and is shown explicitly 
in [DP12]. The bound on the third derivative of Ep tCT in (1.35) and on the first and second 
derivative of if/po- in (1.36) are new and are proven in [DP12]. 

It turns out that the properties stated in Proposition 1.1 are not quite enough for our pur- 
poses. Scattering theory for several electrons requires much more detailed information about 
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the electron's localization in space than scattering of one electron and photons. This informa- 
tion is contained in regularity properties of the momentum wave functions of the vectors i/^o-. 
Let us express if/pa- in terms of its m-particle components in the Fock space: 

to*r = {fp^mm,, (1-39) 

where e Ll ym (R 3m ,d 3m k), i.e., each fp a is a square-integrable function symmetric in m 
variables from K 3 . Let us introduce the following auxiliary functions: 

C(k u . . . ,k m ) := f| cAXbF ^ W , jc^a-eo)" 1 *, 0<eo<l, (1-40) 
\J \k\ i/z 

where c is some positive constant independent of m, cr, P and A within the restrictions specified 
above. Finally, we introduce the notation 

J{ run :={ke^\n<\k\<r 2 ), (1.41) 

where < r\ <r%. Now we are ready to state the required properties of the functions : 

Theorem 1.2. Fix < a < 1/2 and set P max = 1/6. Then there exists Aq > s.t. for all 
PeS = S Pmax , A e (0, A ] there holds: 

(a) Le? {fp^meNo be the m-particle components of tj/p^ and let ttj^ be defined by ( 1.41 ). 
Then, for any PeS, the function is supported in mJ" k . 

(b) The function 

S x 3 (P; fci,..., fc w ) ^ , . . . , £ m ) (1 .42) 

is fw/ce continuously differentiable and extends by continuity, together with its deriva- 

— xm 

tives, to the set S x J^ ,*,. 

(c) For any multiindex fS, < < 2, the function (1.42) satisfies 

< -J=l*/l" W ^(*i,...,*«), (1-43) 
Vm! 

|3j#U*i,...,*i»)l < -^(^-f'^i,...,^), (1-44) 

IVfyJpU*i. •.•.*»)! < ]— * fer^i,-..,^), (1-45) 

where the function Aq h-» /I? /zo^ f/je properties specified in Proposition 1.1. 

Parts (a), (b) of Theorem 1.2 and estimate (1.43) in (c) can be extracted from [Fr73, Fr74, 
Fr] or proven using the methods from these papers. The key new input are the bounds (1.44), 
(1.45) on the first and second derivative w.r.t. P with their mild dependence on the infrared 
cut-off cr. These bounds require major refinements of the iterative multiscale analysis from 
[Pi03] and they constitute the main technical result of the companion paper [DP 12]. 
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1.3 Scattering theory 

In this section we outline the construction of two-electron scattering states along the lines of 
Haag-Ruelle scattering theory [Ha58, Ru62], following some ideas from [Fr, A173]. Since 
we are interested here in the infrared-regular situation, we set 1/2 > a > 0. Such a will be 
kept fixed in the remaining part of the paper. (The infrared-singular case a = is much more 
involved due to the infraparticle problem and will be studied elsewhere). 

Let Cq(S) be the class of twice continuously differentiable functions with compact support 
contained in 5. Following [Fr, A173], for any h € C?(5) we define the renormalized creation 
operator 

oo 

ft a {h) := V — d 3 pd 3m kh(p)f^(k)a*(k) m rj*(p - k). (1.46) 
±d Vm! J 



m=0 

In this expression we use the following short-hand notation, which will appear frequently 
below: 

f^(k):=f^(k u ...,k m ), (1.47) 

a*(k) m := a*{k l )...a*{k m ), (1.48) 

k := h + ■■■ + k m . (1.49) 

It is shown in Lemma A.2 that rf^Qi) and /^(/ji)?)^^), for h\,ii2 € Cq(S), are well defined 
operators on C. (Since fja-(h) := (if^Qi))* is obviously well defined on C, we obtain that ff^Qi) 
is closable). 

Now let O e and Of be the vacuum vectors of r(t) e ) and r(bf), respectively. Then Q. := 
Q e ® Qf is the physical vacuum in 7Y. As shown in Lemma 2.3 below, 



:= = n* J d'PKP)^, (1.50) 

that is ifjh,a- is an element of the subspace of renormalized single-electron states "Hi,^, defined 
in (1.38). Since we assumed that 1/2 > a > 0, we obtain from Proposition 1.1 (c) that there 
exists the limit 

f h := lim ?7*(fc)Q - n* d 3 Ph(P)t// P . (1.51) 

cr^O J 

Clearly, iph belongs to the renormalized single-electron space "Ki of the Hamiltonian H. 

Let us now proceed to the construction of two-electron scattering states. We fix some 
parameter yo > 4, which will be fixed in our investigation, choose some y € (4, 70] and 
introduce a time-dependent cut-off 

cr t := K/t r (1.52) 
for t > max{l,/c-}. Next, we choose h\,h2 e Cq(S), with disjoint supports and set 

h ut (p) := e- iE "%(p), ie{l,2}- (1-53) 
Now we are ready to define the two-electron scattering states approximants: 

%MM := e ifft C ( (/*u)C ( (M". (1-54) 

We will show that the limit of ^tMM ex i sts as t — > 00 and can be interpreted as a physical 
state describing two independent excitations. This is the content of our main result concerning 
scattering theory, stated below. 
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Theorem 1.3. Fix 1 /2 > a > 0, y > 4. Let A e (0, Aq], where A' Q > £j sufficiently small. 
Then, for hi, lt2 € Cq(5) w/?/z disjoint supports, the following statements hold: 

(a) Le? o~ t - K/t 7 , where y € (4, yo]. f/jere exists the limit 

Km : = }!Z e iH % t (hi,t)€ t (h 2 ,t)n (1-55) 

and it is called the two-electron scattering state. It is independent of the parameter y 
within the above restrictions. 

(b) Let h2 , V F^, h , be two scattering states. Their scalar product has the form 

(KM'K' v h'^> = ^h^h\){^h 2 ,^h' 2 ) + {ip h{ ^h'Mh 2 ,^h\)- (1-56) 

Proof. Part (a) follows from Theorem 2.1. Assumptions (2.2), (2.3) of this theorem are verified 
in Proposition 3.4. Assumption 2.4 follows from Propositions 2.2, 3.1, 3.6 and Corollary 3.5. 
Part (b) of the theorem follows from Proposition 3.4 and Proposition 1.1 (c).D 

We remark that an essential ingredient of the proof of convergence in (1.55) is the disjointness 
of velocity supports of the functions h\,h2, defined as 

V{hi):={VE p \pesuwhi}, ze{l,2}. (1.57) 

This property follows from the assumed disjointness of the supports of h\,h2 and the invert- 
ibility of the relation S 3 p ^ E p , guaranteed by Proposition 1.1 (a). 

We note that the vector i/o,, given by (1.51), is determined uniquely by the function h e 
Cq(5). Moreover, vectors of the form tA/,, ®. s </%, where h\,h2 e C^(S) have disjoint supports, 
and ® s is the symmetric tensor product, span a dense subspace in *Ki ® s 'Hi . Thus we can 
define the wave operator W + : 'Hi 8^ 'Hi — > 'H as follows 

1^0^®,^):=^. (1-58) 
By Theorem 1.3 (b), this map is an isometry. 

Standing assumptions and conventions: 

1. The parameters P max = 1/6 and 1/2 > a > are kept fixed in the remaining part of the 
paper. 

2. The maximal coupling constant Aq > is as specified in Proposition 1.1 and Theorem 1.2 
for the values of P max and 1 /2 > a > fixed above. It remains unchanged in Sections 3 
and 4. It is readjusted to a possibly smaller value A' Q > in the last step of the analysis 
in Theorem 2.1. A' may depend of yo but not on y e (4, yo]- 

3. #;i > ^0 will denote positive functions of Aq e (0,Aq], which may differ 
from line to line, and have the property 

lim = 0, lim 5\ = 0. (1.59) 

Such functions control the infrared behaviour of our estimates listed in the statement of 
Theorem 2.1. By reducing the coupling constant from Aq to A' we make this behaviour 
sufficiently mild. 
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4. We denote by y € (4, yo] the parameter which controls the time dependence of the (fast) 
infrared cut-off i.e., cr, = K/t 7 . The parameter yo is kept fixed in the remaining part 
of the paper. This parameter appears also in the definition of the slow infrared cut-off 
cr s = /c(o7/c) 1 ^ 8ro) in the proofs of Lemmas 4.1 and 4.2. 

5. We will denote by c, c' , c" numerical constants which may depend on S , Aq, sq, k, a, yo 
and functions h\,ii2 but not on cr, t or the electron and photon momenta. The values of 
these constants may change from line to line. 

6. We will denote by (p, q) i-» D(p, q), (p, q) \-> D'(p, q) smooth, compactly supported 
functions on R 3 x R 3 , which may depend of S ,Aq, eq, k, a but not on cr, t or photon 
momenta. 

7. We will denote by k - (k\, . . ., k m ) e R 3m a collection of photon variables. A lower or 
upper index m of a function indicates that it is a symmetric function of (k\, . . . , k m ). For 
example: 

f m (k):=f m (k u ...,k m ). (1.60) 

Similarly, we set a*(k) m :- a*{k\) . . . a*(k m ). We note that the order in which the compo- 
nents of k are listed is irrelevant, since they enter always into symmetric expressions. 

8. We separate the electron variable pel 3 and the photon variables k € R 3m by a semi- 
colon. For example: 

G m (p;k) := G m (p;ki,...,k m ). (1.61) 

9. Two collections of electron and photon variables p e M?,k e R 3m and q e R 3 ,r e R 3 " 
are separated by a bar. For example 

F m>n (p; k | q; r) := F m>n (p; k\,. . . ,k m \ q; r\, . . . ,r n ). (1.62) 
10. Given k - (k\, . . . , k m ) we write k :- k\ + ■ • ■ + k m . 



2 Main ingredients of the proof 

Theorem 2.1 below gives the existence of scattering states based on some assumptions proven 
in the later part of this paper: Property (2.2) follows from Proposition 3.4. To establish prop- 
erty (2.4) we derive a formula for time-derivatives of scattering state approximants in Propo- 
sition 2.2 in the later part of this section. The section concludes with Lemma 2.3 which gives 
a direct integral representation of the vectors tph,a- = f)*o-(h)Q- We used this result already in 
(1.46) above. 

In our proof of convergence of the scattering state approximants (1.54) we will vary time 
t and the infrared cut-off cr independently. Let us therefore fix h\,h2 e C^(S) with disjoint 
supports and introduce an auxiliary two-parameter sequence 

% tCr := e m %(h u WAh2, t )n. (2.1) 
Now we are ready to state and prove the main result of this section. 
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Theorem 2.1. Let h\,hi £ Cq(S) have disjoint supports, let *F fi(7 . be given by (2.1). Suppose 
that for A e (0, Aq], infrared cut-offs cr, o~' s.t. o~ < o~' < k and t > max{l, k] 

(%,o-',%,cr) = (^h u a-',^h u a-)(^h 2 ,a-',^h 2 ,a-) + R(t,CT,0 J ), (2.2) 

* c Mic^ +i(rTa \ {23) 

C /rr"^ 4 ^ - 1 1 \ / 1 1 \ 

mA * o^(— + " + 7o^) + c<t1 ' Sa ° + C(T T + 7^)) (2 - 4) 

Then one can choose A' e (0, Aq] s. t. for A e (0, A' Q ] there exists the limit 

Km = }^^„ (2.5) 

where o~ t = K/t 7 . This limit is independent of the choice ofy e (4, yo]. 

Proof. We assume that ti > t\ > 1 are sufficiently large so that cr t2 < o~ tx < 1 . We write 

Wn,a t2 - Vtism 11 * l|VF ' 2 .-2 - + ll^i^ - ^ II- ( 2 - 6 ) 

Concerning the first term on the r.h.s. of (2.6), we note that the bound in (2.4) implies 

|ft^|| S -y2 + 5 -^) + «rf*, (2.7) 
for some e > 0, depending on yo, but independent of io- Now we estimate 

-'2 

it! ' 1 ( r 



WVt^-VtutrJ < ^ dt\\d^ t , at2 \\<CO-^\og{t 2 lh) + cl g J myo) +CCr^t 2 2 



t2 

A/3 



1 1 l/2-<$. 1(1 2 1 ^2 

^ l/Oro) fc 2 ^/ fi 

.i/j ,1/3 g// 

* c (ir + (ir) )' (2 - 8) 



< co£+c 

J/3 .1/3, 



where in the third step we made use of the fact that t~ a logfeAi) is uniformly bounded in 
t% > t\ > 1 for any a > and chose /to sufficiently small (depending on yo) to ensure that 
e' > and that 8^ + l/(4yo) < l/(3yo). In the fourth step we chose Aq again sufficiently small 
and exploited the fact that y > 4 to ensure that 1 > e" > 0. 
As for the second term on the r.h.s. of (2.6), we get 

ll^l,tr (2 _ ^1,07, II 2 
= (™, ,o-, 2 , >cr(2 ) + (Tfj >0 - (] , ,a- tl ) - 2Re( v P, 1 >(T>2 , *P f , >0 - (] > 
= (*Ph l ,a-, 2 ,^h l , < T, 2 ){^h 2 , C r, 1 ,^h 2 , l T, 2 ) + R{ti,0- t2 ,0- tl ) 

+(tffh u o-,, , ^hua-n )(^h 2 ,<T,. , fih 2 ,<T t . > + R(tl,Cr tl ,0- tl ) 



-{fthuo- t ,,fth l ,a t Xfih 2 ,o- tl >&h 2 ,o-, 2 ) + R(h,o- h ,o- h ). (2.9) 
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Thus making use of Proposition 1.1 (c) and (2.3), we obtain 

Wh^-VturJ 2 Cof + \R{h,(T t2 ,(T t2 )\ + \R{h,(T t{ ,(T tl )\+2\R{t { ,(T t2 ,<T h )\ 



1 i\ 1 



U ( 2 '2 



Thus, by choosing Aq sufficiently small, we can find < Si , 62, S < 1 s.t. 

/I tf 1 if 



where in the second step we used that ti>t\. Consequently 



11*^-^11*4*- ). (2.12) 



In view of (2.8), we get 



4 ftl'\e, 



ll^-^rjl<c£ (-*-)'' (2.13) 



for < Si < 1. Let us now set *P(f) := and proceed as in the proof of Theorem 3.1 of 

n-1 



[Pi05]: Suppose fi <t%< t n , +l . Then we can write 



\mt 2 ) - v(t x )\\ < ( Yj h^ +1 ) - + 11^) - 

4 - / 1 \ 4 1 

- C EZ( ( «3-l/3))- C Z'l ,7 1 _(l/ ?1 )2e,/3- (2 - 14) 

/ — 1 /C— 1 J I — 1 

Since the last expression tends to zero as t\ — > 00, we obtain convergence of t \-> *P(0- 

Finally, let us show that the limit *F+ fe is independent of the choice of the parameter 
y e (4, y ]. Let 4 <y' <y and let cr, - x/t 7 , cr' t = K/t 7 ' so that cr t < cr' r We will show that 

lim||1V, -%^\\ = 0. (2.15) 

t—. >oo ' 

Similarly as in the first part of the proof, it follows from formula (2.2) and from Proposi- 
tion 1.1 (c) that 

Um(%^,%^) = Um(%^,%^,) = lim^;,^;) = (i/r hl ,ilf hl )(i/f h2 ,i(r h2 ). (2.16) 
This concludes the proof of (2.15). □ 
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In Proposition 2.2 below we derive a formula for dt¥ tt0 -, appearing in assumption (2.4). For 
the purpose of this derivation we introduce a suitable domain: First, we fix / e No, r\, r 2 > 
and consider vectors of the form 

% un = V — d 3l pd 3m kF lM {p;k)Tf{p) l a*(k) m £l, (2.17) 

where F; >w € Ll ym (B*[ x S^ n ,d 31 pd 3m k) i.e., are square-integrable functions, symmetric 
(independently) in their electron and photon variables and supported in each electron (resp. 
photon) variable in a ball af radius r\ (resp. r 2 ). Moreover, the norms of F/ m satisfy the bound 

m 

||F,, m || 2 < — (2.18) 
Vm! 

for some c > 0, independent of m, which guarantees that the vector (2.17) is well defined. Now 
we set 



2):=Span{Tp^|Z6N ,ri,75>0}, (2.19) 

where Span means finite linear combinations. This domain is dense and it contains C. We 
show in Lemma A.2 that H^, H e , H f , H* /c and fj*(h), h € Cg(5), are well-defined on D and 
leave this domain invariant. 

Proposition 2.2. Le? h\ , h 2 6 Cq(S ) /jave disjoint supports and let ^ t ,a- be given by (2.1 ). Then 
there holds 

+ iffc{h u W a Q%Fl} + {1 ~ 2}. (2.20) 

The operators and are defined on C by 

H* := J d 3 pd 3 kv%(k)7i*(p + k)a(k)7i(p), (2.21) 
%r ■■= J d 3 pd 3 kv<L{k)rf{p-k)a*{k)r 1 {p), (2.22) 

where vZ(k) := A^jf and hf(p) = (E p ^ - E p )hi(p). 
Proof. We compute 

d t % >a = e^iHrf^h^flUHdO. + e iH %(d t h lt t)fia(h 2tt )n + e^ih^d^n 

= e mt iHC(h u )ftUh2,t)^ ~ e iH %(h2, t WMEh) t )n - e iH %{h lt )fi*AKEh 2 ) t )n, 

(2.23) 

where (Ehj)(p) := E p hj(p), i = 1,2. The first term on the r.h.s. above is well defined by 
Lemma A.2. The equality 

(<Wm, ? ))C(M" = ?AdM,,WAh2,t)n (2.24) 
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can easily be justified with the help of Lemmas B.4 and 3.3. Now we note the following 
identity, which is meaningful due to Lemma A.2: 

iHfj^h^fjUh^a = i[[H, fjUhi, t )l vAhjm + fj*AhuWvAh2, t )n 

As for the first term on the r.h.s. of (2.25), we note that [H&, ff a (hi,t)\ and [H?,fj^(hi,t)] are 
sums of products of creation operators and therefore commute with %.Qi2,t)- Thus we get 

i[[H, f,l{h u )l €(fo, t m = fA)]. €(h2, t m. (2.26) 
As for the second term on the r.h.s. of (2.25), we obtain 

+fjUihuWA(Eh 2 ) t )a. (2.27) 

Here in the first step we applied Lemma 2.3 and in the last step we made use of the fact that 
the operator 

Kr-= fd 3 pd 3 kv'L(k)r ] *(p + kMk)r 1 (p) (2.28) 

annihilates ^(/i2 >f )il due to the fact that vZ is supported below the infrared cut-off. As the last 
term on the r.h.s. of (2.25) can be treated analogously, this concludes the proof. □ 



Lemma 2.3. Let h € Cq(S) and let us consider the vector i/^ j0 - := f)l.(h)Q.. Then 

i/fh^ = Tl*r dPh{P)ifj P ^. (2.29) 
Consequently, H^t/f^ - ifr Ea -h,a-> where {E a h){p) := E pa h{p). 

Proof. The m-particle components of if/f,,o-, with the electron's variables in the configuration 
space representation, have the form 

^(x, k lt . , . , k m ) = J d 3 p e^hip + k)f™ +ka {h ,..-,k m ) 

= ^m e ' ik - X \ d 2 pe i P x h(p)f™ a {h,...,k m ). (2.30) 
Consequently, the m-particle components of n(i^/ 10 -) are 

n(ifr h , <r ) m (P,k 1 ,...k m ) = ^ J d 3 xe- iPx f (Ppt? px h(p)j% r Qc 1 ,...,km) 

= h(P)f™ a (k u ...,k m ), (2.31) 

which concludes the proof. □ 



16 



3 Vacuum expectation values of renormalized creation 
operators 

In this section we estimate the norms of the terms appearing on the r.h.s. of (2.20) in order to 
verify assumption (2.4) in Theorem 2.1. We also derive estimates on the norms of scattering 
states, required in this theorem. A crucial input is provided by the non-stationary phase analysis 
in Section 4, which in turn relies on the spectral information from Theorem 1.2, proven in 
[DP 12]. 

In this section we will use the following definitions: 

G Um {q; k) := e'^hiiq)/^), i e {1,2}, (3. 1) 

where hi , h 2 € C^(S ). (We stress that G,> are f-dependent, although this is not reflected by the 
notation). Moreover, we set 

B* m (Gi,m) := J d 3 qd 3m k G i<m (q; k)a*{k) m rf{q - k). (3.2) 

Another convention, which we will use in this and the next section, concerns contractions 
of creation and annihilation operators. It is explained in full detail in Lemma B.2, so it is 
enough to illustrate it here by a simple example. Let us set h = 3, fh - 2, n = 2, m = 3 and 
consider photon variables r e M? n , k e R 3m , r e R. 3 ", k e R 3m . Then the expectation value 

<Q, a(rfa(kfa*(r) n a*(k) m a) 

= (Q, {a(h)a(f 2 )a(r 3 )} {a{k y )a{k 2 )} {a*(n)a(r 2 )\ {a* {ki)a* (k 2 )a* (k 3 )}Q) (3.3) 

is a sum of (m + n)\ terms resulting from all the possible contraction patterns. Let us consider 
one of them: 



(O, {a(h)a(f 2 )a(r 3 )} {a(h)a{k 2 )} {a* ( n )a* (r 2 )} {a*{k l )a*{k 2 )a*{h)}Q). (3.4) 
i I i I 



Given this contraction pattern, we make a decomposition r = (r, r), where r = (?\,r 2 ) are the 
r-variables which are contracted with some r-variables and f = f 3 is the r-variable contracted 
with a ^-variable. In the case of the ^-variables we have k = k, since both k\ and k 2 are 
contracted with ^-variables. In this situation we say that k is empty. Similarly, we have r - r 
with r empty and k - (k, k) with k = (k 2 ,k 3 ) and k - {k\}. 

3.1 Double commutator 

It turns out that the behaviour of the double commutator on the r.h.s. of (2.20) is governed by 
the decay of the functions 

F n!,f 2 i^ r\p\k) := (n + 1) J d 3 r n+ i Va(r n+ i)G hn+ i(q + r n+y ; r, r n+ i)G 2jn {p - r n+ i ; k), (3.5) 

where q, p e R 3 , r e R 3 ", k e R 3m . For any such function we define an auxiliary operator 

B* n , m {Fn!,f 2 )-= j d 3 qd 3 p J d 3n rd 3m kF^q;r\p;k)a\ryW(k)"Y(p-k)q\q-^. (3.6) 

Key properties of expectation values of such operators are given in Lemma B.3 below. Now 
we state and prove the estimate on the double commutator which gives rise to the first term on 
the r.h.s. of (2.4). 
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Proposition 3.1. There holds the bound 

/cr ff/(4ro) i \ 

wvmi, c(h u )i w*ym * ^(— — + * + ^ T7 ^ y ). 0.7) 

Proof. We compute 

([[HI fil{h u )l fcih^tW, [[flf, ff a (h u )l ft a {h 1>t m) 

J] . = ([[H^,B* R (G lfi )],Bt(G2,mm, [[HlB*(G hn )],B* m (G 2 , m m), 

m,n,m,n<EN ^™\n\m\h\ 
m+n=ih+h 



A / -1/1/1 1 r 1 1 j/M 1 11 1 



m,n,m,neN ^m\n\m\h\ 
m+n=m+n 

(S;-i, m (^££) + S;_ 1 ,„(F^))Q) ! (3.8) 

where in the last step we made use of Lemma 3.2 and the operators B* n n ( ■ ) are defined in (3.6). 

Given families of functions Gi un ,Gi 2tn ,Gj un ,Gj 2 , n of the form (3.1), with i\,i 2 ,j\,j 2 6 U>2}, 
we define 

C{Gi x ,Gi 2 ;G j x ,G j^) 

:= ^ / „ ^-.^^ iU^ (3-9) 
ym'.(n + J iimifw -<- i >' 



m,n,m 



yfmMn + l)!m!(« + 1)! 



Then we obtain from (3.8), taking into account that B* } m (F^f 2 )H = B* , „(F^ 2 f >)0 = 0, 

||[[/^C(^)]X(^)]^I| 2 -C(Gi,G 2 ;Gi,G 2 ) + C(G 2 , d; G 2 , d) 

+ 2Re(C(G 1 ,G 2 ;G 2 ,G 1 )). (3.10) 

In view of definition (3.1), C(G 2 ,Gi; G 2 ,Gi) can be obtained from C(Gi,G 2 ;Gi,G 2 ) by a 
substitution (h\,h 2 ) — > {h 2 ,h\). Thus it suffices to consider C(G'j, G' 2 ; G\, G 2 ), where 

G'^te;*) := e- iE "%(q)f^(k), i e {1,2} (3.11) 

and (/zp/ij) € {(/*i,/i 2 ), (/i 2 ,/ii)}. We recall from Lemma B. 3 that 

P^S m+n 

x ffif 2 (P -k + r,r,k\q + k-W)+ l?f\q + k-r,r,k\p-k + r,k r)).(3.12) 



The notation £, fc, f, f is explained in Lemma B.2 and at the beginning of this section. Now we 
obtain from Lemma 4. 1 

i /cr a/(4yo) 1 \ 1 

l^te r I p; ft)| < + crt + —^-J—DCp, q )CWM, (3.13) 

where (p, g) h-> £)(p, g) is a smooth, compactly supported function. From this bound we get 

\F^-(q; r \ p; k)F^f\p - k + r,f,k\q + k- r,k,r)\ 

1 in-"l {Ayo) 1 \ 2 1 „ 



< 



(T 
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and similarly 

n,m V 

1 /n- S/ ( 4 >' o) i \ 2 1 „ 



|F^ G %; r | p; k)F°'f 2 (q + k-r,r,k\p-k + r,k,r)\ 



(—— +crt+ D'(p,q) ^_^ zg^) 2 g» 2 , (3.15) 



where D' is again a smooth, compactly supported function. Making use of (3.9),(3.12) and the 
last two bounds we get 

\C(G' v G' 2 ;G l ,G 2 )\ 

1 /cr^ 4 ^) 1 \ 2 y (m + n)! 

o-^ol t °~ fi(r l /(4ro)) Zj y m i(„ + \)\fh\{n + 1)! 

x — 1= ^8% 
ym'.n'.m'.n'. 

l fo sK4vo) i \ 2 /K t \ 4A2c2 C io-" 1 ^ 1 \ 2 



where we made use of Lemma 3.3 and definition (1.40). This concludes the proof. □ 

Lemma 3.2. For any n,meNo there holds 

[[Hf, fl£(Gi„)], B* m (G 2 ^m = B:_ hm (F^)Q + B* m _ ln (F^ n )Q, (3.17) 

where F^g is defined in (3.5) and we set B\ m (F G _\f n 2 )Q. = B\ n (F G _\^)Q. = 0. 
Proof. First we compute the inner commutator on C: 

We note that 

/i 

[a(w)n*(u + w)j](u), a*(r) n rf(q - r)] - ^ 6(w - ri)a*(r u ) n ~ l n*(u + w)rj{u)rf{q - r) 

i=l 

+a(w)a*(r) n 6(u -q + r)rf{u + w), (3. 19) 

where a*(r u f^ 1 =a*(r x )... a'fa-iKfa+l) • • • a*(r n ) for n > 1 and a* fa.)" -1 = for n = 0. 
Since Gi >n is symmetric in the photon variables, the contributions of the first and the second 
term on the r.h.s. of (3.19) are given by 

[H[, B* n (Gi n )]i 

:=nf d 3 qd 3(n - l) rd 3 r n d 3 uG Un (q; r, r„)v^(r„)a*(r) n_ y (h + r n )n(u)j]* (q - r - r„){3.20) 

and 
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respectively. Now let us compute the first contribution to the double commutator: 

[WlB* n {G u y\uB* m {G 2 , m m 

X [a\rf- l rf{u + r n )q{u)rf (q - r - r n \ a\k) m rf{p - k)]Q 

x a*{rf- l a*{k) m r 1 *{p - k + r n )rf(q - r - r„)Q. (3.22) 
By changing variables p — > p - r n and q — > q + r n , we get 

[[Hi B;(Gi >n )] i , B* m {G 2 , m m = Bl_ l m (F^)Q. (3.23) 
The second contribution to the double commutator has the form: 

[[Hl,B* n {G hn y\ 2 ,B* m {G 2 , n m 

= J d 3 qd 3n rd 3 w J d 3 pd 3m kG Un {q;r)Va(w)G 2 , m (p;k) 

X [a{w)a*{rfrf\q -r + w), a(k) m q*(p - k)]Q 
= m J"d 3 qd 3 p J d 3 "rd 3(m ~ l) kd 3 k m Gi A {q;r)Va(k m )G 2im (p;k,k m ) 
x a\rfa{k) m - l ri*{q ~r + k m )rf(j> -k- k m )Q. (3.24) 
By changing variables q — » q - k m and p — > /? + we obtain 

[[fl* B:(G 1>n )] 2 , B;(G 2 , m )]« = B* m _ u (F^)Q, (3.25) 
which concludes the proof. □ 

Lemma 3.3. There hold the estimates 



2 =«ioM^r, 

m,n,m,heNo 



m+n=m+n 



m,n,m,neNo 
m+n=m+h 



m+n=m+n 



where g™, k* are defined in ( 1.40) and we set by convention g^ = 0. 
Proof. By definition of the functions g™ 



\\C\\\ < {Ac) lm {\og{Kjcr)) m . (3.29) 
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Thus we get 



^— ' „ m\n\m\n\ V / — L mini I mini 

m,neNo 
m+n=m+n 

^--L \cr) 



m,n,m,neNo m,neNo m,»eNo 

m+n=m+H rh+n=m+n 



Let us now prove (3.27): 



S £( 2 ^)^ <1 „ g( . / .,r« 

^— ' \ ^— ; ' m!«! / m!«! 

m,«eNo m,«eNo 

m+fi=Hi+n+l 

^— m\n\ \cr) 



m,n,ih,n£tto m,n,m,heNo 
m+n=m+n rn+n+l=fh+n 



The bound (3.28) is proven analogously. This concludes the proof. □ 



3.2 Clustering of scalar products 

In this subsection we consider clustering properties of scalar products of scattering states ap- 
proximants. We will study the expression 

(% <a ,,%,<r) = <fl, f) a ,(h' 2t )fi^(h\ t )f}UhiMAht)^)- (3-32) 
Recalling from (1.46) that the renormalized creation operators have the form 

oo 

ff a (h) = Y— d 3 p d 3m k h(p)f^(k)a*(k) m rf(p - k), (3.33) 

we obtain that (3.32) is a sum of integrals over vacuum expectation values of the electron and 
photon creation and annihilation operators. Contractions of the electron operators give rise to 
two terms 

, 1=1 1 

(X^^t,^ := <«, ^WWi^W, (3-34) 

i F=i I 

<1^, 1V><*> := <fl, ^(^P^'^'iPC^i^Cfe^n), (3.35) 

which we call 'direct' and 'exchange', respectively. We emphasize that the contractions in 
(3.34) and (3.35) above do not involve the photon creation and annihilation operators. The 
contractions of photon variables are the subject of the remaining part of this discussion. 

As for the direct term, we distinguish two types of photon contraction patterns. The first 
type has the form 

i 1=1 1 

<n, a(k) m a(f) n a*(r) n a\k) m Q.), (3.36) 
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where the monomials of creation and annihilation operators above come from the respective 
renormalized creation operators in (3.34). By (3.36) we mean that all the ^-variables are con- 
tracted with ^-variables and all the r-variables are contracted with r-variables, and that m = m, 
n - n. (That is, k and f are empty in the terminology introduced below equation (3.4)). 
Clearly, for fixed m, n there are mini such contraction patterns and, as we will show in the 
proof of Proposition 3.4, after summation over m, n they give rise to the first term on the r.h.s. 
of (3.38), which is (tyh! v o J , x lfh\,a)k x l>w v o J - l x lfhx,a)- Contraction patterns of the second type are 
those for which k or f are non-empty. It will be shown with the help of the non-stationary 
phase analysis from Lemma 4.2 that the resulting terms contribute to the rest term R(t, cr, cr') 
on the r.h.s. of (3.38), which eventually tends to zero in the proof of Theorem 2.1. 

As for the exchange term (3.35), we consider again two types of photon contraction pat- 
terns. Contractions of the first type have the form 



i 1 1 I 

<Q, a(kfa(rfa*(r) n a*(k) m Q.), (3.37) 

i.e., r and k are empty. These contraction patterns give rise to the second term on the r.h.s. 
(3.38), which is (^h[,a J ^h2,cr)( , ^K l ,a J ^h 1 ,(r)- Contraction patterns of the second type are those 
for which r or k are non-empty. The non-stationary phase analysis from Lemma 4.3 gives that 
these terms contribute to R(t, cr, cr'). 

After this overview, we are ready to state and prove the main result of this subsection which 
yields estimate (2.3). 

Proposition 3.4. Let k > cr' > cr > and let ^q-, *FJ a be two scattering states approximants 
given by (2.1) with momentum wave functions h\,hi and h' v h' 2 , respectively. (We recall that 
supp hi n supp h.2 = and supp h' x n supp h' 2 = 0. However, supp hi n supp h'., i,je {1, 2} may 
be non-empty). Then 

(%,o->^t,<r) = (^h' v CT',^h U (T)(^h' v <T',>Ph 2 ,a) 

+(^h' l ,<T',^h2,<r)(^h' 2 ,a-',^h l ,o-) + R(t,CT,0-'), (3.38) 

where the rest term satisfies 

\R(t,cr,o-')\ < C(h,h')-^(^—^— + (cr') 2 «). (3.39) 

Here C(h,ti) := cpillip2lll Zj8,j3 2 ;O<|0il+|/82l<i ll^'^'ilUII'^ 32 ^!! 00 an d me sum extends over 
multiindices f3\, fc- 

Proof. Let us set 

Gi, m (q;k) := 

for i e {1,2}. Now we can write 

(Y t(T „%^= V — j_ — (H, B h (G\ s )g ft (G^ ft )g;(Gi,„)g;(G 2>w )Q). (3.42) 

m,n,m,nm ^™lnlmlhl 
m+h=m+n 



e ~iE q tU 



(3.40) 
(3.41) 
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Making use of Lemma B.4, we obtain 

(Q.,Bn{G\ fi )Bf h {G'2^)B* n {Gi^ n )B* n {G2,m)^) 



= 2 f J~ d 3n rd 3m kG: J! (q:r)G 2jl Ap:k) 

P^S rn+n 



xyGijjiq + k-f;r, k)G 2 m(p -k + r,k,f) + G lR (p - k + r,f, k)G 2fh (q + k-r;k, r)j, (3.43) 

where the notation in (3.43) is explained in Lemma B.2. 

Let us denote the summands involving the first term in the bracket on the r.h.s. of (3.43) 
by / — Let / ~ ~ be such summands coming from permutations for which k or f are 

J m,n,m,n m,n,m,n b * 

non-empty (cf. the discussion below formula (3.4)). We note that there are (m + n)\ - m\n\ 
such permutations. Let (^p.,,^^) be the contribution to (3.42) involving all the summands 
n mn- Making use of Lemmas 4.2 and 3.3, we get 

1/11 W • ^J 2 " 2 

K V ,T„)<», <_ C(M -,_L(I_ T ±_ + (0 -)(^ 
Clearly, this expression contributes to R(t, cr, cr'). 

a/1 \ { 1 \ y 

Now let r - . denote the expressions r - - coming from permutations for which k and f 
are empty. (We note that there are m\n\ such permutations and that in this case m = m,n = h). 
We obtain from Lemma B.l that 

Crt* = f d^qd^rG^iq^G'^ r) J d 3 pd 3m kG 2 , m (p;k)G 2 ^(p;k) 

1 1 

n! m! 

Let (*P' ,, v F, fr ) (1) be the contribution to ( V F' ,, % (r ) involving all such /*■ ^ Since the sum 

x f,cr" *> cr ' x t,cr" l ' cr ' o m,n,m,n 

over permutations in (3.43) gives the compensating factor m\n\, we obtain 



= —(a,B n {G\ h )B* n {G hn )Q.)(a, B m (G' 2gl )B* m (G2^)n). (3.45) 



O^T^ - V - F z^=<«,B«(G , li j)B*(G 1>n )Q)<fi,B fi (G^)B;(G 2 , m )«) 
nvMM ymlnlmlnl 

- {^h\,o-'^h u a-){^hf 2 ,a',^h2,a), (3.46) 

where in the last step we compared definition (3.2) of B*(G !>n ) with definition (1.46) of the 
renormalized creation operator and recalled that ^. >cr = if^Qi^Q.. 

The analysis of the summands involving the second term in bracket on the r.h.s. of (3.43) 
is analogous, so we can be brief: Let us denote such summands by / ~ - and let I ~ . 

«=> ' J m,n,m,n m,n,m,n 

be the summands coming from permutations for which k or r are non-empty. We denote by 
( x P^ cr ,, x P f tr > (2) the contribution to (3.42) involving all ^ nfhh - By Lemmas 4.3 and 3.3, we get 



K^,*V> (2) I < Cih,h')^r{- t -^ + (3.47) 



1 H 1 

cr A o 

This part contributes to the rest term. 
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*(2) (2) A 

Now let v - - be the expressions r nfhh coming from permutations for which k and f 

are empty. Denoting the corresponding contribution to (3.42) by ( x P' t ^, x I't,<r) \ we get from 
Lemma B.l that 

(X^,^^ = {rlfh'^^a-W^a'^huo)- (3-48) 
This concludes the proof. □ 

The above proposition has a simple corollary which gives the last term on the r.h.s. of esti- 
mate (2.4). 

Corollary 3.5. Let h\,h2 € Cq(S) have disjoint supports and let h^iq) - (Eq <cr — E q )h2(q). 
Then 

||C(/*v)€(^)"H * co- 1/2 ^(l + ^-^). (3.49) 

Proof. Making use of Proposition 3.4, we get 

llfjUhMA^ltMl 2 = (tfrkur, tyhxMh^, + R(t, <r, o-), (3.50) 

where R(t, cr, cr) satisfies 

W,cr,o-)\ < C (h,h')^(]-±^ + (o-^), (3.51) 
C(h,h) = cllAiHiHAJIl! J] ll^fcillooll^^lloo. (3.52) 

ySl,yS2;0<|/?i|+|/?2l<l 

Making use of Proposition 1.1, we obtain that C{h, h) < e'er and therefore 

\R(t,o-,o-)\<cc-^(l + ]-^). (3.53) 

Now making use of Lemma 2.3, we get ||^fc liCr || 2 = \\h\\^ and IIiAa^o-II 2 = PP^- Exploiting 
Proposition 1.1 again, we get ||tA)£,o-ll 2 < ccr 2 . Consequently, due to the constraint on the 
supports of hi and hi, we obtain 

\\fl*Ahi,t)mi t )m 2 < cV^l + 7^^)). (3-54) 
This concludes the proof. □ 



3.3 Contribution involving 

& I,cr 

This subsection is devoted to the term involving H? on the r.h.s. of (2.20). The following 
elementary proposition, which relies on Lemma B.5, gives the second term on the r.h.s. of 
(2.4). 

Proposition 3.6. Let H? be defined as in (2.22). Then there holds the bound 

< co-^o . (3.55) 
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Proof. We rewrite the expression from the statement of the proposition as follows: 

<C(^i,^ I c , (T C(^, t )n,C(^i, t )^i, ( .C(%,f)"> 

= y 1 z (n,B R (G 2t n)(HlJ*B m (G 1 , A )Bl(G hn )Hl^BUG2,m)^ (3-56) 

m,n,ihMN Vm!n!m!*! 
m+n=m+n 

Now Lemma B.5 gives 

(Q, B ri (G2,fi){Hl (T )* Bj h {G\jn)B* n {Gi^ n )H[ a B* m {G2,m)Q) 

= Yj f d ^ d3 P \ d 3n rd 3m kG u (q;r)G 2 , n (p;k) 

xl I d p v-(p) G 2 ,n(p -p-k + r,r, k)G Uih (p + q + k- r,k,r) 

+ P~Hp2,«(<7 + k-r,r, k)G h dp -k + rX ?)). (3.57) 
Making use of Theorem 1.2, and of definition (3.1), we obtain the bounds 

|Gi,„(<7; r)G 2 , m (p; k)G 2fi (p - p-k + r,r, k)G]_^(p + q + k-r,k,r)\ 

1 D{p,q)C{r) 2 g™{k)\ (3.58) 



^m\n\m\h\ 

\G hn (q; r)G 2 ,,„(p; k)G 2fi (q + k-r;r, k)G ltifl (p -k + r;k, r)| 

< \ D(p,q)&(r) 2 C(k) 2 , (3-59) 
ym\n\m\h\ 

where (p, q) \-> D(p, q) is some smooth compactly supported function independent of a, t. 
Consequently, the r.h.s. of (3.57) can be estimated by 

c\\v£\l(m + n)l - L= \\g n jl\\g™\\l (3-60) 
ym\n\m\h\ 

Substituting this bound to (3.56) and making use of Lemma 3.3, we get 

\\f)Uhu)HiJUh2,,W\ 2 < c\%\^f ' • (3-6D 
Exploiting the fact that llv^lb ^ ccr, we conclude the proof. □ 



4 Non-stationary phase arguments 

In this section we derive non-stationary phase estimates which entered into our analysis in 
Subsections 3.1 and 3.2. The spectral information from Proposition 1.1 and Theorem 1.2 is 
crucial for this part of our investigation. 

Lemma 4.1. Let Gj <m , i e {1,2}, be as specified in (3.1) and let F^J^ 1 be defined as in (3.5) 
i.e., it has the form 

Fnf 2 ^ r\p;k) = (n+\) J d*rvz(r)e- i{E ^ +E r- f)t h x (p - f)h 2 (q + f)^^, r)f™_ ?>(r (*),(4.1) 
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where h\,h-i € Cg(5 ) have disjoint supports. There holds the bound 

1 /rr s /( 4 ')'o) 1 \ 1 

where D e C~(R 3 x R 3 ). 

Proof. Let us introduce the slow cut-off cr s := k(o-/k) '^ 7o \ which clearly satisfies cr < cr s < k. 
Let^ € C°°(R 3 ), < x ^ 1> be supported in S\ (the unit ball) and be equal to one on < B\- S for 
some small < e < 1. We set^i(^) = x(k/o~s),X2(k) = 1 - XiO 1 ) an d define 

r\Plk):= \ ^rv^Xj^e-'^^'hip - f)h 2 {q + f) 

We set V(q,p, f) :- E q+f + E p -~ r and note that by disjointness of the velocity supports of h\, 
h 2 , the condition h\(p - f)h2(q + f) ± 0, together with Proposition 1.1(a), implies that 

\V f V(q,p,r)\ >s' >0 (4.4) 

for some fixed e' . Thus we can write 

{-it)\V~ r V(q,p,r)\ 2 

Now we define the function 

\VfV(q,P,r)\ 2 

where x K £ C£°(R- 3 ) is equal to one on S K and vanishes outside of a slightly larger set. We note 
that, by Proposition 1.1 (a) for any multiindex (3 s.t. < |/?| < 1 

\eP f J{q,p,r)\<D{q,p\ (4.7) 

where (q, p) i-> D(g, /?) is a smooth, compactly supported function. Moreover, for < < 2 

where X3 is a smooth, compactly supported function, independent of cr. In addition, for < 
< 2 we obtain from Theorem 1.2 (c) 



< 1-1--^^). (4.10) 



Now using the Gauss Law we obtain from (4.3) 



, (4.11) 

r=crn 
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where n is the normal vector to the unit sphere and dQ(n) is the spherical measure (that is 

n = (sin 9 cos tp, sin 9 sin tp, cos 9) (4.12) 

and d£l(n) = sin 9d9d<p in spherical coordinates (9, <p)). Let us consider the first term on the 
r.h.s. of (4.11). Let/i be the corresponding integrand. Making use of (4.7)-(4.10), we obtain 



1 cD(q,p)xi(f)\r\ 2rt 



-£(/■)«*(*), (4.13) 



where in the last step above we made use of the fact that \r\ < cr s in the region of integration. 
Now let li be the integrand in the boundary integral on the r.h.s. of (4.11). Making use, again, 
of bounds (4.7)-(4. 10), we get 



\h\ ^ — f== — ~ ~ — 8<r( r )8 A k > 



1 cD{q,p) X 3(m 2 ° 
V^W cr 6 * \r\ 2 
1 cD(q,p) X3 (crn)cr 2 " 



Vm!«! cr^o cr 2 

Thus (4.11), (4.13), (4.14) give 



(4.14) 



\F^( q ;r\ P ;k)\ < ' CD{ l P)ia f\ lir) g "Xk), (4.15) 
ymlnl cr"* t 

which is the first contribution to the bound (4.2). 

Now we consider the contribution above the slow cut-off. In this region we will have to 
differentiate by parts twice, so a direct application of the non- stationary phase method would 
result in third derivatives of the function V(q,p,f) := E q+ ? + E p -~ r . However, our spectral 
results do not include the existence of third derivatives of p i-> E p , but rather the bound 



tfpEpjA < — (4.16) 
cr * 

for any multiindex /3 s.t. < |/?| < 3 (cf. Proposition 1.1 (a)). To be able to exploit this bound, 
we introduce an auxiliary function Va-(q, p,f) := E q+ f^ + E p -- r ^ and rewrite F~ u 2 {q; r \ p; k), 
defined in (4.3), as follows 

F 2,n G m 2 (<?; r\p;k) = J d'rv^XiWe-^^h^P - f)h 2 {q + r) ■ f^Jr, f)f^ h(T {k) 

+ R(crt). (4.17) 

The rest term R(tcr) satisfies 
\R(tcr)\ < J d 3 rva(r)x2(m ~ e^^-^^hip - r)h 2 {q + r)f^r, r)f™ f (T (k)\ 

< c -55-£(*)£(r), (4.18) 
cr * 

where we made use of the fact that \V(q, p, f) - Vo-(q, p, r))\ < ccr (Proposition 1.1, (a)) and of 
the bounds (4.8), (4.9), (4.10). This gives the second contribution to the bound (4.2). 
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Let us denote by ^V'bjh tne ^ rst term on tne r - n - s - °^ ( 4 -17)- We wn ^ estimate this term 
with the help of the method of non-stationary phase. Similarly as in the first part of the proof 
we note that by disjointness of velocity supports of hi, hi, the condition h\{p - f)h2(q + f) + 
implies that 

\V?VAq,P,r)\>s" >0 (4.19) 

for some fixed s" , independent of cr € (0, k]. (Here we made use of the fact that 5 3 p i-> E p ^ 
is strictly convex, uniformly in cr. Cf. Proposition 1.1 (a)). Thus we can write 



(-it)\V T -VAq,pJ)\ 2 



e -iVAg,P,f)t = " ~,n ■ ( 4 - 2 °) 

We define the function 



JAq,P, r) := ^f^'^L hip ~ ~r)h 2 {q + f) X \f), (4.21) 

which is analogous to the function / introduced in (4.6) above. Making use of (4.16) and of 
(4.19), we get for all multiindices f3 s.t. < |/?| < 2 



\4Mq,P,~r)\<^-. (4-22) 
By integrating twice by parts in the defining expression for F. 1,c?2 , we get 

xV ' ' ( Jirr^4 V? ' P> ~ r >^X2{f)f n q tL^ ?)f p % «,(*))). (4.23) 

\\VrVcr(q,p,r)\ l ' ' I 

In view of (4.16), (4.19), the function 

, ~ V- r Vo-(q,p,r) 

\yfV<r(q,p,r)r 

is bounded by c/<t 6a °, together with its first derivatives, on the support of (q,p, f) i-> h\{p - 
r)h 2 (q + r)x K (f). Thus we obtain from the bounds (4.8), (4.9), (4.10) and (4.22) that the inte- 
grand / in (4.23) satisfies 

' 1 " *' P) *^-M>a». (4-25) 



V^W(°"s) 2 o*o |r| 2 



where in the second step we made use of the fact that (l-e)cr s < \r\ in the region of integration. 
Thus we get from (4.23) that 



(a-r\vk)\ < 1 cD{q ' p) 1 

ym\n\ cr ' (o s ) z 



\F^Jq;r\p;k)\ < — ( 4 -26) 



which gives the third contribution to (4.2). The factor (n + 1), appearing in (4.1), can be 
estimated by 2" and incorporated into the constant appearing in the definition of □ 
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/(I) 

m,n,m,h 



Lemma 4.2. Let f^J ^^be defined as follows 

4!U« := f d ^ d3p \ d 3 "rd 3m kG Ut (q;r)G 2 , m (p;k) 

X G lR (q + k-r;r, k)G 2fh (p -k + r;k, r), (4.27) 

where Gi t „, G' 2 „ appeared in (3.40), (3.41), the notation k - (k,k), r - (r, r) is explained in 
Lemma B.2 and we consider a permutation in (3.43) for which k or f are non-empty. Then 
there holds 

CmI * (ifi + Zmfi) + C^) 25 ) ^L^ OteT'lg + \\g n X)(\\8TX + ll^li)X4.28) 

where C(h,h!) :- cpi||i||/i 2 |li ZjSiA^lSil+l&l^i II^^HcoH^ 2 ^ 00 awd ram extends over 
multiindices /3\, f3 2 . We note that for k (resp. f) non-empty we have m + (rasp, n + 0) and 
there always holds m + n = m + h. We set by convention g^ 1 = 0. 

Proof. By inserting the definitions of G J>n , G' ~, i € {1,2}, we obtain 

X*l(P - * + D^(<7 + * " Df2Ar)fp,Ak)f q+ l_-i,Ar> k)f p _ k+b(r ,(k, r)-(4.29) 

Since the expression on the r.h.s. of (4.29) and the bound (4.28) are invariant under the substi- 
tutions k <r* r, p q,h\ <-> /j 2 , /j'j <-» /j 2 , m <-> n, m <-» n it suffices to consider the case of 
non-empty L Similarly as in the proof of Lemma 4. 1, we introduce a slow infrared cut-off. Let 
us set cr's :- k(o-' I k) ii ^°\ which clearly satisfies o J < cr' s < k. Let^ e C°°(R 3 ), < # < 1, 
be supported in Si (the unit ball) and be equal to one on Si_ s for some < s < 1. We set 
X\(k\) := x(k\/o J s ), ^2(^1) := 1 _ ^fi(^i)» where &i is the first component of k, and write for 
7 e{l,2} 

m,n,m,n I I 

Y.h x (p)h 2 (q)h x (p -k + r)h 2 (q + k- f)f^(r) f^(k) 

v — /"? v ——in w 

x(*/*i)/V*-£^ *)) fp-l + i,A k > f y ( 4 - 3 °) 

Let us first consider (4.30) with j = 1. We conclude from Theorem 1.2 and the definition of 
the functions g£ that 

1 ^iA)^(^Wg":(r,%"(l,r) 



1 

^m\n\m\h\ I&1I 3 



1 trfiWk[«(ii) 2 liif" „, , 2 



Vm!n!m!n! |&i| 3 



j&WsT'tfT, (4.31) 
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where we decomposed k - Qi\,k'),k - ik\,k') and in the last step we made use of the fact that 
g%(k,r) < g™{k,r) and g^,\r,k') < g^ l {f,k') for <x' > cr. Substituting (4.31) to (4.30), we 
get 

ym'.n'.m'.n'. 

Let us now consider (4.30) for j = 2. We set w := k?_ — r and define 

V(h,p,q,w) := E_ h+p _ w + E h+q+w . (4.33) 

We note that, by disjointness of velocity supports of h\, h 2 , the condition h\{-k\ +p-w)h 2 (ki + 
q + w) + implies that 

\V h V(k u p,q,w)\ >s >0, (4.34) 

for some s' independent of k\,p,q and w. Thus we can write the following identity 

,„r s Vjt V(ki,p,q,w) • Vr e w ^ - p * w)? 

it\V h V(k u p,q,w)\ 2 

Now we define the function 

V~ k V(k u p,q,w) _, v _, v 

J{k\,p,q,w) := — — /liC— A:! + p - w)h 2 {k\ + q + w). (4.36) 

|V i[ V(&i,/?,g,w)| 2 

We note that, by Proposition 1.1 (a), for any multiindex /? s.t. < |/?| < 1 there hold the bounds 

|fl£/(*i,.p,?,w)| < c (h\,h' 2 ), \8^jiii)\ < (4.37) 
where co(/z' 1 , ftj) has the form 

c (h\,h' 2 ) = c J] ll^'^IUII^IU, (4.38) 

/?i,/32;0<IAI+|/3 2 |<i 

where /?i,/?2 are multiindices. Moreover, we obtain from Theorem 1.2 (c) 

|/^(r)| < -^=C{r), (4.39) 
Vn! 



s Mzkf^- f) - (440) 



m!V(cr') ^ 

l#C*2(*i)7^_^fr%^(*))l * n={7^ + ^T^k)g7(k). (4-41) 
fcl Vn!m! Ho" ) ^ °V 

Now coming back to formula (4.30) and integrating by parts we obtain 

?(1)(2) 



, p, 9 , *) fpAWp-hlAK r)\ (4.42) 
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Making use of the bounds (4.37), (4.39), (4.40), (4.41) and of the fact that gfc(f, k) < g%{r, k), 
since cr' > cr, we estimate 

*m™ 1/1 1 \ll^lllllfellico(^,/io) , o 

* Vt^U V jlIlOl- (4-43) 

? V(cr') ov ym\n\m\n\ 

Exploiting (4.43), (4.32) and the fact that cr^ = k(o-' /k) 1/( ^°\ we get (4.28). Here we estimated 
trivially 

ii^nliigr 1 ^ + rtiiszil ^ + ik&iixiisr'iia + wcwl) (4-44) 

to obtain an expression in (4.28) which is symmetric under the substitution m <-> n. □ 
Lemma 4.3. Le? 7 ;n n ^ . fee defined as follows 

X G un (p - k + nr, k)G 2M (q + k- fjk, f), (4.45) 



: = J ^qd 3 p j '^"rd""kG lM (q:r)G 2jll (p:k) 



where G^ m , G'. , i e {1,2}, are defined in (3.40), (3.41), the notation k - (k,k), r = (f, r) is 
explained in Lemma B.2 and we consider a permutation in (3.43) for which k or f are non- 
empty. Then there holds 

* ()(4r + TTT^)) + (^)-^=(\K-% + rtG^lg + ll^li),(4.46) 

where C(h,h') :- cpi||ip2||i 2/?! ,/? 2 ; o<L0i |+|/3 2 |< i lld^ilUII^ 2 ^ 00 an ^ sum ext ends over 
multiindices f}\, We note that for k (resp. f) non-empty we have m + (resp. n + 0) and 
there always holds m + n — fh + h. We set by convention g~ l = 0. 

Proof. By inserting the definitions of G,\ m , G' we obtain 

Z,fnfi= f f d 3n rd 3m ke iiE ^ E ^ e ^ P -k + r E ^% 
xh\(q + k- r)h 2 (p -k + ^(r)/^^ f)f^(k)fl_ k+t(r ,(r, k). (4.47) 



/(2) 



We note that by substitutions h\ <r* h' 2 ,m <r* h,k <r*k,r <r* f we obtain formula (4.29). Now 
the statement follows from Lemma 4.2 and the fact that (4.28) does not change under the above 
substitutions. □ 



A Domain questions 



Lemma A.l. There exist constants < a < 1 and b > s.t. for any Y e C (n) there holds the 
bound 

Wf^W^aW^W+bm, (A.l) 
where H^ 1 ' = Hi\ C („), H£ = Hf r \ C („) and the constant b may depend on n. 
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Proof. Let us use the form of the interaction Hamiltonian appearing in formula (1.25). We 
have = + , where 

n „ 

H *,(n) I d 3 kVa(k)e ikXi a(k) (A.2) 

and H c 1 ' {n) = (H\ M )* . Let us set C n (k) := E" = i ^ and compute for some ¥ e C (,i) 
II^W^I < J rf 3 iv5(A:)||C„(*)a(A:)T|| < n J d 3 *v 5 (*)||a(A;)1'|| 

< «|h- 1/2 v^|| 2 (^,// f »P)2 <i||flfrT||+n 2 ||a>- 1/2 vs|||imi, (A.3) 

where in the last step we anticipate that (A.l) should hold with < a < 1. 

Let us now consider the creation part of Hj . Making use of the canonical commutation 
relations, we get 

wh^w 2 = J d%d% v ^k 1 )^k2)(c n (k 1 ra\k 1 )%c n (k2ra(k 2 m 

< \\H^ n) n 2 + n 2 \\v«\\ 2 2 \m\ 2 , (A.4) 
which, together with (A.3), concludes the proof. □ 

Lemma A.2. The domain D, defined in (2.19), is contained in the domains of H, H e , Hf, H^ c , H^, 
and ff(h), h e C^(S). Moreover, these operators leave D invariant. 

Proof. Let Yj 1 '' 2 be a vector of the form (2.17). Then 

= Y^L-L fd 3 p'd 3l pd 3m 'k'd 3m kF[ m ,(p';k')F hn (p;k) 

x n*(p')r,*(p) l a*{k') m 'a(k) m Q. 

- V -— rf 3(/+1) prf 3 *ikF A (p;ikV(p) /+1 fl*(jk)*n, (A.5) 
*=o v^! J 



where F[ m ,{p', k') := h(p' + ^)f^ +k Jkf), m:=m + m',k:= (k, k'), p := (p, p') and 



Fm(P; ~k) = Y, ■= ' ^ {F\ (ni _ m) F lm ) sym {p- k), (A.6) 

where the symmetrization is performed in the p and k variables separately. By Theorem 1.2, 
Fj , satisfies the bound (2.18), and therefore 

c'" 

H^alb < ~r=, (A.7) 
Vm! 



for some constant c. Hence ffeihyS 1 , 1 ' 2 is well defined and belongs to £> 
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Next, we note that 



oo 

HJl?<* = V _ d 3 'pd 3m k + • • • + Oipi))Fi^p\ k)rf{p) l a\k) m a, (A.8) 

00 1 f 

= y d 3l pd 3m k {oj{k x ) + ■■■ + co(k m ))FUp; k)if( P ) l a*(kTn. (A.9) 

f-i Vm! J 



Due to the support properties of Fi >m (p; k) these vectors are well defined and belong to D. 

Finally, we consider the 
stricted to *H® has the form 



Finally, we consider the operators H^ c . We recall that the interaction Hamiltonian re- 



H ? := 2j f d3kv a(k)(e ikx, a(k) + e~ ikx, a*(k)) (A. 10) 

(=1 ^ 

and we express in terms of its m-photon components, i.e., 

{% ur2 f' m \ P l, ...,p,;k u ...,k m ) = F hn ( Pl ,...,pr,ki,..., k m ), (A. 11) 

Now we can write 

(Hl% i ' r *f> m >( Pu ..., P i-,k 1 ,...,k m ) 

= VmTT f ^v^)2](^ 1 ''" 2 ) ( '- m+1) (^ 1 ,..., A -fc,..., /7/ ;^,^ 1 ,...,^ ! ). (A. 12) 

It is easy to see that for some constant c, independent of m, 

||(^^^)(t«)|| 2 < -£=. (A. 13) 



/m! 

Similarly, we obtain that 

(^ 1 ' r2 ) ttm) (p 1 ,...,^^i,...,^) 

= — V V v^W^f'^iPu ■ ..,Pj + k, . . ., P r,h,. . .,k u , . . .,k m ), (A.14) 

where fe^ means omission of the i-th variable. This gives, again, a bound of the form (A. 13). 
Since the case of is analogous, this concludes the proof. □ 



B Fock space combinatorics 

In Lemmas B.l, B.3 and B.5 below we verify the identities first for G, jm ,G'. m ,F nfn ,F' nm of 
Schwartz class and then extend them to square integrable functions using Theorem X.44 of 
[RS2] 

Lemma B.l. Let G m , G' m e L 2 (R 3 X R 3m ) be symmetric in their photon variables, see ( 1.61 ). 
Let us define as operators on C 

B*JG m ):= [ d 3 pd 3m kG m {p;k)i 1 *{p-k)a\k) m (B.l) 
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and B m (G m ) := {B* m (G ,„))*. Then there holds the identity 

(a,B m {G' m W m {G m )Q) = m\ J d 3 pd 3m kG m (p;k)G m (p;k). (B.2) 

Proof. We compute 

(n,B m (G' m )B* m (G m )Q) 
= J d 3 pd 3m k J d 3 p d 3m k! G m (p;k)G' m (p';k')(Q,a(k') m j](p' -V)rf(p-k)a(k) m Q) 

= J d 3 pd 3m k J d 3 p' d 3m k' G m {p;k)G' m {p';k')6{p - k- p' + V)(Q,a{k') m a*{k) m Q.) 

= J d 3 pd 3m k J d 3m k' G m {p;k)G' m {p - k + V;k')(£l,a{k') m a*{k) m £l) 

/„ m 
d 3 p d 3m k I d 3m k' G m (p; k)G m (p - k + k'; k') ^ \\ s <- k P® ~ k 'i) 
J peS m i=l 

- «^W"*<WH^;*), CM) 

where S m is the set of all permutations of an m-element set and in the last step we exploited 
the fact that G' m is symmetric in its photon variables. □ 

Lemma B.2. Let n,m,h,m e No be s.t. n + m = h + in. Let us choose 

r = (n,...,r n )eR 3n , k = (k u ...,k m )eR 3m , (B.4) 

f = (n,...,f s )eR 3 ", * = (*!,...,**) eR 3 * (B.5) 

and define the sets 

C :={l,...,n}, C; := {n+ \,...,n + m], (B.6) 

C»:={l,...,n}, Ci :={n+l,...,n + m}. (B.7) 

(TVote ?/io r ? is the complement ofC n in {1, . . . ,n + m\. Similarly for C'~). Let S m+n be the set 
of permutations of an m + n element set. For any p e S m+n we introduce the following notation: 

r '■= (ri)(i, P (i))ec„xc-„, r '■= ('"/)(!, P (i))eC„xq> (B.8) 
^ ■= (k~n\i,p(i))eC' n xC^ k :- (ki- n \ itP (i)) e c' n xCn> ( B -9) 

so that r = (f, r), k = (k, k). Similarly, 

r := (r p (i))(i,p{i))zc n xc~ n , r := (? p (,))( 1 >(;))€QxC ii , (B.10) 
k := {kp(i)~h)(i, P (i))eC' n xC^ k := (&p(i)-n)(i,p(i))eC„xCS (B.ll) 

50 ?/jo r ? f = (r, f) a«<i £ - (fc, (7f{ / 1 (i,p(i)) e C n x C„) = % then we say that f « empty, and 
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analogously for other collections of photon variables introduced above). Finally, we define 

S{r-r) := f] ^ ~ f K®>> (B " 12) 

(i,p(i))eC n xCa 

S(r-k) := Y] 5( ri -~k p(i) - n ), (B.13) 

(i,p(i))eC„xq 

:= Y] Kh-n-fpif)), (B.14) 

(/,p(0)6C^XCs 

5{k-k) := Y] m-n-~kp(i)-*)- (B.15) 

(;,p(0)ec;,xq 

r/jen ?/iere /zo/cfa 

(Q, a(rfa(k) ih a*(r) n a*(k) m n) = ^ s (?~ *Wf ~ ~ r)S(k - k). (B.16) 

P^S m+n 

Proof. Let (vi v n+m ) - (n, . . . , r n , k\, . . . , k m ) and (vi, . . . , v n+m ) = (h h,h, ■ ■ ■ , km)- 

There holds 

m+n 

<Q, a(rfa(kfa*(r) n a*(k) m £l) = £ ft 6(v J ~ 9 */>) 

psS m +„ j=\ 

= ^ <5(ri - v p (i)) . . . <5(r„ - v p ( n) )(5(fci - v p( „+i)) . . . 6(k m - v p(n+m) ) (B.17) 

= z ( n ^ -*>«))( n 

pe5,„ + „ ((>(i))eC„xCH (/,p(j))eC„xq 

x( f] %_ n -r p(0 )V P] <%■_„- £ p(0 _ n )l (B.18) 

(i>(0)eC£xC s (i,p(i))eCixq 

which concludes the proof. □ 

Lemma B.3. Let F n>m , F' n m e L 2 ((R 3 xR 3 ") x (R 3 xR 3m )) be symmetric in the photon variables. 
Let us introduce the following operators on C 

K,m( F n,m) := J d 3 qd 3 p J d 3 "rd 3m k F n , m (q;r\ p;k)a*(rf a*{k) m rf(p - kjr]*(q - r) (B.19) 
and set B n ^ m (F n j n ) := (B* nm (F '„,„,))* . There holds 

(B n<m (F' m yn,B ntm (F n , m yn)= £ J d 3 qd 3 p J d 3n rd 3m kF n , n (q;r\p;k) 

X\F' hm (p -k + r,f,k\q + k-r,k,f)+ ~F- nih (q + k-r;r,k\p-k + r,k,f)^ (B.20) 

for n + m = h + m. Otherwise the expression on the l.h.s. is zero. Here S m+n is the set of 
permutations of an m + n element set and the notation k, k, f, f is explained in Lemma B.2. 
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Proof. We compute the expectation value 

(Bl^F'^B^iFn^Q.) (B.21) 

= J d 3 qd 3 pd 3 qd 3 p J d 3h rd 3m M 3n rd 3m k% m {q\r\p^ 

x(8(q - p +k~ r) + 5{q - q - f + r))8{p + q - p - q) 

x<0, a{ffa{kfa{r) n a{k) m Q). (B.22) 
The last factor is non-zero only if n + m - n + m. Then 

<fl, a{ffa{ky h a*{rfa*{k) m a) = J] 8{r - f)S{f - k)5{k - f)8{k - k), (B.23) 

where we made use of Lemma B.2. Thus the r.h.s. of (B.22) is a sum over p € S m+n of terms 
of the form: 

J d 3 qd 3 pd 3 qd 3 p J d 3n rd 3m k F '^(q; f,k\p;k, r)F n<m (q; r\p;k) 

X 6(p + q - p - q)(S(q - p + jc- r) + 5{q - q -k + r)) 
= J d 3 qd 3 p J d 3n rd 3m kF n , m (q-,r\p\k)(F' m {p-k + ntk\q + k 

+ F'hfiiil + l-nr,k\p-k + r;k,r) ) J (B.24) 

which concludes the proof. □ 

Lemma B.4. Let Gi >m , Gj , G2, m > Gj m 6 L (R 3 x R 3m ) be symmetric in the photon variables. 
We define, as an operator on C, 

B* m (Gi, m ):= j ' d 3 pd 3m kG Um {p-k)rf\p-K}a\k) m (B.25) 

and set B m (G^ m ) = (B* m (Gi^ m ))* . There holds the identity 
(Q, Bjj(G^ fi )Bft(G 2 ^)B*(Gi )n )B^(G2, m )Q) 

= f d 3 qd 3 p J d 3n rd 3m kG Un (q;r)G 2 , m {P> k ) 

G l R {p - k + r;r, k)G 2 ^(q + k- r,k,f) + G i h {q + k-r;r, k)G lih {p - k + r;k,r) I, (B.26) 

for n + m - h + in. Otherwise the expression on the l.h.s. is zero. Here S m+n is the set of 
permutations of an m + n element set and the notation k, k, f, r is explained in Lemma B.2. 

Proof. Follows immediately from Lemma B.3. □ 

Lemma B.5. Let G hm ,G 2 , m e L 2 (R 3 x R 3m ) be supported in l 3 x(iel 3 | > cr} xm and 
symmetric in their photon variables. There holds the identity 



(Q, Bn(G2fi)(Hl a )* Bm(G\,m)Bl(G\, n )Hi a B* m (G2,m)Q) 



psS,„ 



x( I d"pv^(prG2,n(p-p-k + r;f,k)G lt0l (p + q + k-r,k,f) 

+ \K\\ 2 2 G~2fi(q + l-r,f, kjG 1Jh (p -k + r,k, F>) (B.27) 
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for m + n - m + h, otherwise the l.h.s. is zero. Here S m+n is the set of permutations of an m + n 
element set and the notation k, k, f, r is explained in Lemma B.2. 

Proof. We compute the expectation value 

= J d 3 ud 3 wd 3 ud 3 w J d 3 qd 3 pd 3 qd 3 p J d 3n fd m kd 3n rd 3m k 

v^(w)v~(w)G 2 ,n(q', r)G\,fh(p; k)Gi A (q; r)G 2 , m (p; k) 

x(fl, i](p - k)rf(u)n(u - w)n(q - r)rf(q - r)rf(u - w)r]{u)rf{p - k)Q.) 

x<0, a(rfa(w)a(kfa*(r) n a*(w)a*(k) m Q,). (B.28) 

We note that 

(O, J](p - k)q\u)q{u - w)rj{q - r)rf(q - r)n\u - w)n(u)n*(p - k)Q) 

- S(p -k_- u)8{p - k— u)(Q., i~i(u - w)i](q - r)q*(q - r)q*(u - w)Q.) 

- 6(p -k_ - u)8{p — k — u)^6(u - w — q + r)S(q -f-u+w) 

+ 6{u -w - u + w)S(q -f- q + r)J. (B.29) 
Let us consider the contribution to (B.28) of the first term in the bracket in (B.29): 

(Q., B^(G2fi)(Hi a .)* Bm(Gi^m)Bl(Gi^)H^ a .B* m (G2,m)^)l 

:= J d 3 ud 3 wd 3 ud 3 w J d 3 qd 3 pd 3 qd 3 p J d 3n fd 3ih kd 3n rd 3m k 

v^(w)v^(w)G 2 fi(q; f)Gijh(p; k)G\ >n (q; r)G 2 , m (p; k) 

x5(p — k — u)6(p - k - u)8{u -w - q + r)S(q -f-u+w) 

x(n,a(ffa(w)a(kfa*(r) n a*(w)a*(k) m a) 

= J d 3 ~qd 3 pd 3 qd 3 p J d 3n fd 3lil kd 3n rd 3m k 

v^(w t )v^(w t .)G 2 ,n(q; r)G]_ M {p; k)G hn (q; r)G 2 ,m(p; k) 

x<£2, a(ffa(kfa(w t )a* (w» )a* (rfa* (k) m Q), (B .30) 

where in the last step we integrated over u,w,u,w and set w* := p—k—q + f,w* :- p—k—q + r. 
Now we consider the expectation value of the photon creation operators: 

(Q,a(r) n a(kfa(wM(w*)a*(r) n a*(k) m Q.) = 5{w, -w -,)<Q, a(ffa(kfa*(r) n 'a '(£f Q>,(B.31) 

for r, k, w, f, k, w in the supports of the respective functions. (Here we made use of the fact 
that |w*| < cr, whereas |r,-| > cr, \kj\ > cr). Let us now substitute the r.h.s. of (B.31) to (B.30). 
Making use of Lemma B.2, we obtain 

(Q Bn(G2fi)(Hi a .)*B rn (Gi^m)Bl(Gi^)Hl a B* m (G2,m)^-)l 

= f d 3 qd 3 pd 3 qd 3 p J d 3n fd 3ih kd 3n rd 3m k 

P^S m+n 

- 1 - q + d^ip -k-q + D G 2,h(q; r)G]_ M {p; k)Gx, n (q; r)G 2 , m (p; k) 

x5(p + q-p- q)5(r - f)5(f - k)6{% - k)5(k - f). (B.32) 
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By integrating over q, r, k, we obtain 

(Q, B s (G2^)(// I C ^)*Bm(Gi i m)B*(Gi,„)//[ cr B*„(G2,m)fl)l 

= \ d 3 pd 3 qd 3 p J d 3n rd 3m kvZ(p-q-k + r) 



&\2 

I a pa qu p I a ra k v—^p — q — k -t- 

peS m + n 

x G 2 ,n(P + q- p;r, k)G\jn(p; k, r)G hn (q; r)G 2 , m (p; k) 



= J] J d 3 pd 3 qd 3 p J d 3n rd 3m kv*(pf 



P^S m+n 

x G 2 ,a(p - p-k + r,r, k)G\^(p + q + k-r;k, r)G\, n (q; r)G 2 , m (p\ k), (B.33) 

where in the last step we made a change of variables p^p + q + k-r. This gives the first 
term on the r.h.s. of (B.27). 

Let us now consider the contribution of the second term in the bracket on the r.h.s. of 
formula (B.29): 

(£2, BfiiG^nXH^)* Bm{G {,m)B* n (Gi, n )Hi cr B* n (G 2 j n )Q.) 2 

:= J d 3 w J d 3 qd 3 pd 3 qd 3 p J d 3h fd 3rn kd 3n rd 3m k%{p -I- p + k + w)v^(w) 

xG 2Jl {q + r - r; r)Gi tl}l (p; k)G\, n {q; r)G 2 , m (p; k) 

x<Q, a(ffa(p -k-p + k + w)a(kfa(r) n a(w)a*(k) m Q,) 

= [d 3 wf d 3 pd 3 qd 3 p J d 3n rd 3 ™kd 3n rd 3m k%{p - I- p + k + w)v^(w) 

xG 2Jl {q + r - r; f)G]_ fil {p; k)G hn (q; r)G 2 ,,„(p; k)8{p -k-p + k) 
X(Q, a{ffa{kfa*{rfa*{k) m a) 

= \\vZ\\j Yj J d VP J d 3n rd 3m kG hn (q;r)G 2 , m (p;k) 

X G 2fi {q + k-r,r, k)G Xfil {p -k + r,k, r),(B.34) 

where in the first step we integrated over u, w, u and in the last step we made use again of 
Lemma B.2. This gives the second term on the r.h.s. of (B.27) and concludes the proof. □ 
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